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ABSTRACT 


An analytic model is developed of baroclinic currents 
measured hy several investigators under the Arctic Ocean ice 
cover. A time-dependent, Boussinesq, quasi-geostrophic 
vorticity equation is developed to describe the flow. 
Separahle normal mode solutions are found which show sta- 
tionary current patterns gualitatively consistent with 
observed data. An initial value problem is then formulated 
using density perturbations at a newly forming lead as a 


mer~cing function. 
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T. INTRODUCTION 


fe OBJECMIVES OF “RESPAR@GH AND REVIPVRP Or PREy TOUS RESr ARGH 

The objective of this study was to develop an analytic 
model of specific small scale haroclinic current patterns 
which have been observed on a recurring hasis heneath the 
Arctic ice cap. The patterns are termed small in that they 
are on a scale less than that of the general Arctic circula- 
tion. In order to develop the model a theoretical investi- 
gation was made of the baroclinic response of the Arctic 
Ocean to perturbations in its internal density structure at 
a newly forming lead. 

It is worthwhile at this point to briefly discuss pre- 
vious research in the area of Arctic dynamic modeling, and 
to describe the baroclinic currents measured hy several 
investigators, in an effort to indicate the relationship of 
mars study to previous work. 

Investigations of Arctic Ocean dynamics can be classified 
by their scale as ene as hy the mode of motion described 
(either internal, surface, or a combined mode). The largest 
scale studied concerns processes which affect the entire 
@eean, such as long term mean ice drift and water circulation. 
im order to practically "study motion on@this scale, spacing 
between grid points cannot feasibly be less than tens of 
kilometers. This necessitates parameterizing some of the 
smaller scale processes occurring hetween grid points, rather 


PidiecemviInG Lome tiem explicitly. Factors which must he 





considered on this scale are the shape and bathymetry of the 
ocean basin, Coriolis effects, wind and i1ceystresce, ec. 

Several theoretical investigations of large scale Arctic 
Ocean dynamics which closely match observed patterns of water 
circulation and ice drift have been carried out in recent 
years. Among these, Camphell [1965] proposed a large scale 
analytic ice drift model and Galt [1973a] developed a working 
barotropic numerical model of overall surfacé water circula- 
tion. Both these models were concerned with surface motion 
or an equivalent barotropic circulation. A complete model 
ei large scale anternal Arctic circulation has not yet heen 
published, but work along these lines is heing carried out 
at NOAA's Geophysical Fluid Dynamics Lahoratory, Princeton, 
New Jersey [Semtner, 1973]. 

The smallest scale of analysis concerns the short term 
local effects of phenomena such as refreezing leads and ice 
cracking. These effects result in dynamic processes such as 
forced convection due to the salt rejection and heat transfer 
during freezing. Grid size of concern here would he from 
Mess than a cm, a microscale, to a few meters. Coriolis 
effects would not be important at the smallest end of this 
eeale but effects such as diffusion would have to be con- 
sidered. Several investigators have done work on the problem 
of forced convection associated with freezing of sea ice. 

Foster [1972] treated the problem analytically after 
having published several earlier papers on the problem of 


memeeccmeenveoct On in the oCccCan. the anitial state inv has 





study was a homogeneous under-ice layer. He then assumed 
density to be a function of salinity only and solved the 
equations of motion and salt diffusion directly associated 
with a refreezing lead. Using a two-dimensional model Foster 
@erri1ed his solution to the point of finding initial vertical 
velocities and length and depth scales’'for convection cells 
during the initial stages of convection. These were ahout 
imc S@C fOr vertical velocities, 808em for horizontal 
@Benvection cell spacing, and several cms for depth of con- 
vective penetration. He did not analytically consider 
further stages of convection nor did he continue to a larger 
scale in which Coriolis forces would become important. He 
did postulate a higher order of convection cells with depth 
penetration of up to 6 meters and vertical velocities of 
about 1 cm/sec. 

Piascek [1970] considered the general problem of surface 
cooling effects on the ocean. He developed a numerical 
model for convection cells on both a small scale (order of 
ems) which neglected Coriolis effects and a larger scale 
feecer of kms) which considered Coriolis. His model con- 
figuration was a rectangular hox directly below the surface 
with an initially homogeneous structure. He considered 
uniform surface cooling or boundary condition and then 
developed a vorticity equation in terms of a vector poten- 
tial. This was numerically integrated over a 2-dimensional 
grid to give steady state temperature and vorticity patterns 
PRcnhinehs Dox. With some Moditication his model could he 


adapted to describe Arctic under-ice convection. 
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Coachman [1966] considered a mass and energy balance on 
the layer of water immediately below the ice cover. He 
assumed a downward flow direction under the refreezing leads 
Gue to sinking of highly saline water and a compensating 
upward flow spread out over the area between leads, as shown 
in Figure 1. His study used mean values of freezing rate, 
temperatures, lead size and distance-hetween leads to deter- 
mine mean vertical velocity values. His scale of analysis 
was necessarily larger than that of Foster in order to arrive 
at mean velocities. His downward vertical velocity values 
fase) ON the order of .01 to .08 cm/sec. The depth of his 
Bpertace circulation layer was 100 meters. if the Arctic 
Ocean is viewed from a large scale standpoint, the net 
result of the investigations hy Coachman, Foster, and Piascek 
would be to show that there is a small downward vertical 
velocity associated with each refreezing lead due to forced 
convection and that there is a local increase in the density 
of the water directly under a refreezing lead. 

A related small scale problem which has been investigated 
is that of the ice formation rate and rates of salt rejection 
and heat transfer at a refreezing Arctic lead. Schaus [1971] 
developed a numerical model of a refreezing lead which treats 
these problems and which could be used in coupling the air- 
sea interaction problem with ice dynamics and with the under- 
mee Convection problem discussed here. 

A vertical velocity similar tothat caused hy the convec- 
Eve process May also arise from the initial cracking and 


1Ge motion during lead formation. This specific problem 
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Figuse 1. Schematic of Convection Model Proposed 
Dy COodenmam | !9eG). 
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has not been previously studied. In this case the theory is 
as follows, assuming ice shears upon cracking over a period 
of a few days as shown in Figure 2. There is a stress 
imparted to the underlying water at the ice-water boundary, 
which in turn causes a convergence of the Fkman layer. This 
im EUrnN results in an impulsive vertical velocity under the 
newly forming lead. The net result is that, for an entirely 
different reason, a vertical velocity similar to that caused 
by convection 1s associated with shear across newly forming 
leads. Both these lead-associated phenomena would be classed 
as small scale effects. 

By interacting with the existing density and flow 
patterns these velocity and density perturbations may have 
some, as yet undetermined, effect on the large scale internal 
flow dynamics of the Arctic. A proposed theory for this is 
Giscussed and developed in part B of this section. Any 
mesultant flow pattern due Eo such an interaction meee 
then be expected on a scale intermediate in size hetween 
the smallest and largest scales discussed; on the order of 
meters in depth and Welkers Ee Kilometers ay lengeEn.  Aeeheos 
retical investigation of baroclinic internal flow in the 
Arctic on this scale has not previously appeared in the 
literature. It is important to note here thatthis is the 
scale of a number of internal current patterns or events 
measured in the AGeCrNe . 

The observed currents appear to be baroclinic in nature 


elrenougy ansatistactory theory for their origin has not yet 
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Figure 2. Proposed Shear Mechanism 
for Lead Pormarren. 
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been developed. An example of the vertical distribution of 
horizontal velocities measured by Galt [1967] under Ice 
mland T-3 is SHO in Figure 3. The corresponding density 
@ectribution is also shown. Similar results found by 
Coachman and Newton [1973] are shown in Figure 4. It can be 
seen that a peak in the horizontal velocity occurs in the 
pycnocline, with a maximum at ahout 150 meters, indicating 
memaroclinic origin. Figures shows the horizontal distri- 
bution of velocities found by Coachman and Newton at 150 
meters. It is a plot of velocity vectors measured from 

one position on the ice while drifting over a given area. 

As can be seen, there is a periodic change in the vatleactiics 
eeeection with an apparent wavelength of about 30 km. The 
reason for classing this periodicity as spatial rather than 
temporal is that the observed currents were found to he 
essentially in geostrophic balance. For this type feature, 
phase propagation speeds will be small compared to advective 
speeds. Figure 6 shows a vertical cross-section of veloci- 
ties as reported by Coachman and Newton [1973]. 

A need thus ast (us for an investigation of the Arctic 
Ocean internal flow dynamics on the intermediate scale 
described. The immediate purposes of this study are to offer 
em explanation of the observed baroclinic currents. A 
further use of such a study is to couple the small scale 
phenomena to the overall circulation, which in turn may aid 
in the development of a three dimensional Arctic circulation 
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Bee “THEORY DEVELOPERSIN FHES Tittets 

A theory was postulated to describe the dynamic effects 
on the underlying ocean of a newly forming Arctic lead. The 
theory developed was begun with the premise that a small 
vertical velocity perturbation or density increase is created 
during the formation of each lead, caused by either refreez- 
ing or cracking, as descrihed earlier. The scale of analysis 
used was tens of kms in the horizontal and the full depth 
Of water (N-3000m) in the vertical. 

Gace water directly helow the Arctic iceé cover can typi- 
cally be expected to consist of a shallow homogeneous layer 
(2@—-50 m), Under which a stably stratified density distribu- 
tion exists. There is also some horizontal slope to the 
isopycnals in order to support a mean flow. The horizontal 
gradient will be several orders of magnitude less than that 
of the vertical density gradient in the pycnocline. Pertur- 
bations such as the vertical velocity and density changes 
occurring under a lead, which is a very small area on the 
assumed length scale, should distort the isopycnals over a 
relatively short period of time and short distance. This 
would act as an impulsive event on an overall time scale 
measured in days. As the isopycnals attempt to return to 
m@meir stable condition, a transient waveform of some sort 
would then appear and propagate out from the lead. 

The first problem would then be to find the form of the 
waves which can exist in a typical Arctic density structure 


and to then determine if these are consistent with observed 
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current patterns. Once this was done, a second problem would 
be an investigation to determine whether the perturhations 

pera lead -could actually tormce these osce1]laki0ons. oneness] 
tigation of these two problems forms the basis of this thesis. 

An i1ndueatbrommet the possable validity Cf tne wore. at any 
point is the degree of correlation of theoretically developed 
waves or currents with real data. A comparison was therefore 
made with the baroclinic currents measured by Galt [1967] and 
Coachman and Newton [1973] once a theoretical current pattern 
was developed. 

Cia Srep Ine SvostAittaGinG sabe MMeory Was  e@uanVves= 
tigate the type of waves which could exist in the given 
density structure and to compare these with known data. This 
corresponded to treating the Arctic Ocean as a semi-enclosed 
system, subject to certain governing equations of motion and 
boundary conditions, and solving for its preferred modes of 
Sssealation. The actual solution to this problem is the 
subject of Section II of this thesis. The general assumptions 
made in formulating the theory and the rationale behind them 
will be discussed ao 

The equations used follow the development of Galt [1973h] 
Wereneromoased On that of Kuo [1973]. The basic difference 
from Kuo's development is in the formulation of the equation 
of state. The system was assumed to be inviscid, hydrostatic, 
and quasi-geostrophic. The equation of state assumes density 
proportional to salinity only. This is a reasonahle approx- 
imation in the region being investigated as can be seen in 


Migises ~awilen saa plot Of Salinity, temperature, and 
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density profiles taken under T-3 along with current measure- 
ments [Galt, 1967]2 

The Boussinesg approximation was also made, which states 
that density variations are significant only when the density 
is multiplied hy gravity. An additional assumption necessary 
to this development is that the Rossby number is small. This 
allowed the expansion of the dependent variables as power 
series in terms of the Rossby number. Equations consistent 
in powers of Rossby number could then be derived. 

From the above assumptions and the basic equations of 
motion a quasi-geostrophic vorticity equation was developed 
in terms of the perturbation pressure (pressure could he 
written as the sum of a mean and a perturbation quantity). 
This equation was solved by assuming a separahle solution 
consisting of the product of a horizontal plane wave and an 
Micnown vertical structure function. This vertical structure 
function was then solved for numerically subject to appro- 
priate boundary conditions. Results of the numerical solu- 
£i0n are compared with the measured currents in Section III 
of this thesis. 

Investigation of the second problem described (whether 
Emme theoretical oscillations could actually be driven hy 
perturbations found at a lead) was accomplished by using 
pressure conditions under a newly forming lead as an initial 
Sendition for the quasi-geostrophic vorticity equation 
leveloped Ania oeculOnmil mA GeSCLIiption Gf this initial value 


problem forms Section IV of this thesis. 
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Ti. DEV EEOEMEN TT OF ie Wave SOLE TLON 


eA nY TTC DEV ABO MENT OF THE VORTICITY EQUATION 

The initial formulation of the problem was carried out 
using a rectangular coordinate system, with positive values 
@eex tO the East, y to the North, and z upward. The govern- 
ing eguations of motion for the flow can be written in the 


following Boussinesgq, inviscid form: 


ee 
u, + uu, + ae ow Ui fv = . Py (1) 
Viet 0 Tea eta wie 4 ti = — os 2} 
t x Me “ ed 
w, + uw. + vw_ + ww_ = - ae choos _— g (3) 
e x y Z De | Oe 


where u,v, and w are the x,y, and 2 velocity components, 
respectively, f is the vertical component of the Coriolis 
parameter, g is the acceleration of gravity, p is the pres- 
muce, Pp 1s the density and P. is the horizontal space and 
time average of the density, which is a eee ee. ZaeOel vy, 
1s@ Be dition to the equations of motion, the continuity equa- 


mheom Can he written as: 


Py, + up + vp_ + we, + De Ve Ie a (4) 


x a4 


and the equation for conservation of salt, assuming a con- 


servative law with no mixing, is: 


SUS pet Vs, se 0 es) 
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EQUdewonS 1s) Seese cen Dyeme aigeenem se VScalted (and. non 
dimensionalized using essentially the same system as was 
used by Kuo [1973]. The non-dimensional quantities (primed) 


eee defined as follows: 





(x,y) = L(x',y'); (ug = Ulu',v'); (6) 
yee izes w= ee 
1 
t= Pl et f= ££"; 
p= p, +t £,ULp_¢'; o = UP aasior 
s = S.(1+es'); 


where L is a characteristic horizontal scale, D is a charac- 
meutrotic depth scale, U is a characteristic horizontal 
velocity and P. and S. are the static pressure and density 
fields respectively. 

In addition the following non-dimensional coefficients 


Biseomael ined: 


Lo ee Rosshy number; (oe) 
O 
ot 
eS wen oa GatetO.Or the Seale length, ti, Eo 
w 4 EReeRO@csny Cerormatton lenge; 
See he, a measure of the perturbations in 
the thermodynamic quantities; 
1 
6 = — yp 7 &@ Measure Clperescegmoiml tay sot tne 
P Ps Stabrveedensity -Lrele; 
J TS 6 
Oe = = yp 7 Aa Measure OF wehbe Sstabllity “Ot ‘tie 
S Shere sealer ye bela, 


SeOStltUttOn freon (6) nto ud) ana (2) then Gives: 
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' iat tal leant 1 LE gree a 
R(uy, tu upity Bet Be uss) flv “ie (8) 


R(viy 


a ee we) tw ° = =m | | (9) 
x yy, Z ¥ 
mma Substitution from (7) into (4) gives: 


(2) 0 


erent at et Ne Lo! Ln 
- Shek Aa wos + Rw Woo) Oo oat 0 eG 


Mreeviding that the hydrostatic relationship for the steady 


state density and pressure fields is used, as given by: 


item, Substitution Gf (6) into (4) Gives: 


2 0 MRS 9 SS (2 MS pee — ee ee 
S 


mir y, SUbDStLItut1On Lrom (6G) into (5) gives: 





6 
ete Se BUS peor DiS ey ae Rai ey eee i w') = 0 (12) 
mK y Z. u 
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Equations (8) through (12) can then he simplified by consi- 


dering limitations which occur in the real ocean: 


> << 1 used in eq. (10) (13) 
ee Peed in cael) (14) 


ena by assuming a simplified equation of state relating 


@ensity to salinity: 


Sareea eu -Cinenotonale torn or the sperturoation Quantities. 


op" = s!; Sree ene Oe ie (15) 
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Usang equa (lee, theswenscriencionalsequaitions)jjecome 


(omitting primes from here on): 


R(u, -|- uu, + a + Rwu,) = os ee aaa (16) 
R(vy } uv, +- ae +4- Rwv ,) ve EON ae a, ers 
d. + Sl + s = 0 (18) 
uoo+v + 0 ob (19) 
x Y J SZ 
eI 
R(s, +. us, + ee + RWS , + i 7s) ae (20) 


At this point the assumption of a small Rosshy number, R, 
will be used and dependent variables expanded in a power 
Sem~1es £Orm; 1.6€., u, v, 6, and s will take the form: 


Ca R a Ca 
n=0 


where the a represents the general dependent variable and w 


Vee pe written as: 


Ww = W + ae He tee cee ey 


+ RW 5 3 


1 


This is consistent with the scaling given by equation (6). 
In addition, the Coriolis parameter will be scaled ona 


beta-plane: 


f 1 + By OZ 


O 


Il 


where 


By < 0(R) (23) 


Using these approximations plus eq. (15) and the series 


(nN 


ie) 





expansion gives the following equations to zero order in 


Rossby number: 


PN =O (24) 
x 

ae mE (25) 
y 

oA = 358 (26) 

Z 

U ty = 0 (27>) 
O O 
x y 


These show the flow to be geostrophic, hydrostatic, and hori- 
memcally non-divergent to a zeroth order approximation. The 


equations to first order in Rossby number are: 








ao Fouou, a ss Cee. oe BY = —94 (28) 
t Xx my. Xx 

Yo, Ue. eee Byu, = Ore 29); 

x By 

on = or i oa (20) 
Z, 

Cet aac Se). = 0 ei) 
Oe Vy S 
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Ses. Vase = sit Noes 0 (2 ) 

1e X 4 


To eliminate the first order pressure terms from (28) 
and (29) the equations are cross-differentiated and sub- 


i@eacted, resulting in: 


(vue -u_ ) + u_(v_ -u_ ) + v_(v. -u_ ) + (v. -u_ )(u_ tv_ ) 
ee on fe Qe, a ce 7 oe . 2 7 


+ hy +: a er + BY = 0) (33) 


ici oop ctitirvon from ({2/) and (31) “into (33) results in: 
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— 


ih _ 
D_ mes), * ey 18, (34) 


meeaition (32) is rewritten in the following form: 


ie 


os | on cso 0.’ ~ ° Sie 





(wP,) a 


Tye 

migseiteution Of (26) imto (35) (noting that = # eres) ) 
ie 

gives: 


Up 


S 
(wi P.) = We 
0 Zz 


Pie 


Cure ve | (36) 
P a P “a Nh 




















and differentiating this with respect to z gives: 



































oe Pie Bie 
CHE oS akc T UVOlB %6 Vols % 
p ZB ZU Ome 2 7x fe ZZ 
Up up 
S S 
sy 4 6 6 _ GS é %o | (37) 
Zz 0 ZX Z, 0 Z! vy 
firme Can them je substituted directly into (34). In writing 


ifeeet te 1S helpful te derine all et the velocities in terms 
of the pressure, assuming a geostrophic relationship (this 


assumption filters out gravity waves). 


(38) 


MPEG hiecce rte Deoaseen that themilact two’ terms in (37) 
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identically cancel elUhen, substitution of (37) and (38) 


into (34) and dropping all zero subscripts gives: 
Se Wee Se = (39) 
where: 


= St Be 
BO Win IY 0. i 2) (40) 


Equation (39) is the quasi-cgeostrophic vorticity equa- 
mom and (40) as the potential vorticity. The first two 
terms on the right hand side of (40) constitute the relative 
vorticity. The third term is the contribution due to the 
VeGtation in the Coriolis parameter, and the last term repre- 
Gets the vorticity contributed by vertical stretching of the 
water column. Equation (39) shows that time-dependent motion 
will result when the motion is such as to create a change in 
Poremticalivorticity. Porsthe anitial preblem studied here, 
the Coriolis parameter, 8, was assumed to he non-zero and 
the resultant motion investigated. It should be noted that 
as the North pole is approached, 8B approaches zero, but that 
motion could still be driven by bathymetric variations. 


Assuming a separable wave solution of the form: 


oe 3 (z)er (kxtly-st) (41) 


leads to the following Sturm-Liouville problem for the verti- 


cal dependence: 





; 
Q 2 Qk 7 
y h Os ee @ 2 = 0 (A) 
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Substituting £romeq2.a6/) sor 1 sand Di listings the faceuehnat 


dp. | 
the Brunt-Vaisala frequency, N, is given by : 7 and 
S 





converting hack to a dimensional form results in: 


ZL 
O Z eee Z. Bk es 
; lo, = : “FL oh Bx | Z = 0 (43) 
N Z 





B. EQUATION USING ARCTIC DFNSITY PROFILE 
ioe basicC otobempescription 

In this initial formulation zero mean flow was 
assumed, thus the initial density distribution was assumed 
to be a function of depth only. The Arctic under-ice density 
profile was approximated by a homogeneous upper layer from 
the surface to z = -h, below which density increased follow- 
ing an exponential curve. In this lower layer density 
increases such that it assymptotically approaches a constant, 
Par at deeper depths. This distribution is similar to that 
proposed by Reid [1948], in his equatorial model, and is 


given analytically by: 


Dames [ (Oz) 


ey 


I od zs 3a0)) 


0 = Dg a (Pa-Pyde (z = = ine) 


which is illustrated in Figure 8 for various values of c and 
h (referred to hereafter as layer depth (h) and e-fold depth 
(c)). Actual density values taken during the collection of 

velocity data are also shown on Figure 8 for comparison. As 
GCanebe Seen; the best data fit is with an e-Lola depth of 


between 100 and 120 meters and a layer depth of 30 meters. 
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2. )Usesonrctucly Detroit hy ae eno lone en Fquations 
The density distribution given by eq. (44) was sub- 


stituted into eq. (43) to produce a vertical normal mode 
eaecion which would apply to this specific density structure. 
To show this, eq. (43) is first expanded in terms of P. and 


mes derivatives: 


a Ps PsPs Z wy ARR | 
2 | Z, + 22 | ommmmanary 12 + ik” +1 epee 80Re (45) 





g 0 ZZ. Z 
sv (op. ) 


The derivatives of p. are found from eq. (44) as: 


— (0 4~Py,) l/c (2th) 


O — 
Sy Cc 
(46) 
_ an Pn) 1 /e (24th) 
Pes = y) e 
Void ce 


wen can then be substituted into eq. (44) to yield, after 


some rearrangement of terms: 


2(p 4-0.) . ; 
, 7 aan’ te tat re 1|, 
ZZ, Oo Cc cl Zz 
.S 
- —L (p,-0,) i? ee Bi f,1/e(etn)g (47) 
ie, C 2 
O 


This can be simplified to: 


Zi - faet/e(24h) + Bz + pel /c (2th) Z = 0 (48) 
TAG 7. 
where: 
BiG Ons) 
Beale 


a5) 





(49) 


OM lie 


2 2 Bk 


<2 aie pee eae o Bs 
Ane : 5 (04 24) + 1° + 3 
Co 


O S 


The solutionsto eq. (48) will then he a set of eigen- 
functions describing the vertical distribution of the pressure 
meanction, >. It is»noted that this equation only applies to 
the region below the layer depth. In the upper homogeneous 
region horizontal velocities will be constant since there is 
Memscratification. Vertical boundary conditions on the 
system are that w, the vertical velocity, is equal to zero at 
both the surface (ice) and hottom boundaries. In scdieen 
velocities are assumed to be continuous across the internal 
boundary at z= -h. If these houndary conditions are met 
and a solution to eq. (48) exists there will he solutions 
for different discrete values of XX, the eigenvalues for this 
equation. 

A value of either k or o may then be assumed and used 
with the eigenvalues to solve for the corresponding values 
Sumponwor k fOr each mode. in order to actually solve for both 
forma O SOme initial distribution for @ must be known and 
iaeched. This is the initial value problem, based on the 
perturbation at a lead, which is the next phase of the study 
ame 1S Giscussea In Seceion iV. 

The problem being discussed here was formulated with 
the eeu eten Ota long=enarrow Teac alveneds:n a NOrPth=souch 


Givcecer ome ewaiue 3 where x, and Yy, are tein Teale eat cre 


e ! 1 
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and North-South dimensions of the lead, respectively. Based 
on this it was then assumed that no changes occur in the y 
Girection (d/dy = 0) which in turn means that 1 is assumed 
to be zero in eqs. (41) and (43) and that u =-d¢/day = 0. 

The equations thus become two-dimensional. A horizontal 
boundary condition at time zero of v = 0 at the origin, 
which is placed in the center of the lead, is then imposed, 


The horizontal velocity can be written as: 


ee Aye ee! 
x 
Sr: 
a kZ(z) [C, cos (kx~ ot) + C,sin(kx-ot) ] 
mem v{(0,0) = 0, C, must equal zero or: 
Vi = Sean) Sabri eo ) (50) 


This means that ¢ can be written as: 


0 —— 7a COs ko Ec} 


From eq. (36) w can he written as: 
ae 
w= a Z Sin (kx-ot) 


which becomes: 
Lo 

aS ay Z Sin(kx-ot) (S29 
N Z, 


when the terms are dimensionalized and proper substitutions 
made for u and Sy. Assuming a value for k is known, an 


expression for o can then be derived from eq. (49): 


35 


-Bk 
5 = (53) 
LD Gp 
re 
J (P57) 


This is similar to the equation for the frequency of a haro- 


tropic Rossby wave, o, = B/k, with the additional term in 


b 


the denominator being due to the stratification. 
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Cas Ka ee Re OU) oe eA oooh AND COMER LSOM 
WITH ACTUAL DATA 
Poe NUMERICAL SOLUTION OF THE VERTICAL MODE EQUATION 
The vertical mode equation (48) was rewritten uSing a 


mememealized depth varrable, 0 > 2° > -1.0, as: 


_DB(z'+h/D) | DB (z'+h/D) 


Z = 


+B)Z - dx 
Tb Td Z 


Z= 0 (54) 


where A! = D7 and D = depth of water. 


Equation (54) was solved numerically on the Naval Post- 
graduate School IBM/360 computer. The solution was found 
Metng a Milne predictor-corrector method, according to an 
algorithm described by Nielsen [1964]. The depth range was 
divided into 2000 increments and then the Milne method used 
meen the bottom, z' = -1.0, upward to the normalized layer 
depth, z' = - h/D. The Milne method section of the numerical 
solution was then tested by using it to solve a zero order 
Peesel equation which had a known analytic solution. Results 
of this test were satisfactorily checked with a table of 
seggel FUME TONS | 

In the relatively thin homogeneous region above the layer 
depth, the slope of the vertical velocity was, held at its 
layer depth value, as was the value of the vertical structure 
function, 2. The output of the numerical scheme was a set of 
2000 values of Z and A, ieeOneuae ese seme! wo FOr aac; 
ee ihe veremeet= Gisteatlioution Of horizontal velocity is 


PuOperelonal co 4 as Was Shown by eq. (50). From eq. (52) it 


can be seen that the vertical dependence of vertical velocity, 
wa LS proportional, to 2 /N*. The numerical scheme could thus 
be used to give the vertical dependence of hoth horizontal 
and vertical velocities. 

ithe bet lLomebeundany cOondmerenge w — Omater' s=s—1.0jawas 
met by setting Z,(-1.0) = 0. This was valid since w is 
proportional to z/N* and since RiaeeniGno) , although quite 
small, is not identically zero. With the bottom houndary 
Semertion thus fixed, the condition w = 0 at z' = 0 could 
only be met for certain values of '. These values are the 
eigenvalues corresponding to the normal modes of oscillation 
of the ocean with the given vertical density distribution. 
Memomacn tO find the courece eigenvalues an initial A‘ was 
assumed. An iterative process was then begun in which the 
numerical method described earlier was used to arrive at a 
value of w(0) for each new value of A'. The initial value 
chosen for A' was then increased in steps until the value 
of w(0) changed sign. At that point a Regula-Falsi method 


(eq. (55) below) was used to find the next A'. 


het a 


ntl ~ W (55) 


This was continued until the value of w(0) was sufficiently 
eose to zero. The value of A" at that point was then con- 
Sidered an eigenvalue. The process was continued by increas- 
ing the value of A' in steps until the ere eigenvalue was 


fLOuUnGd 3m che Same Manner. 
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me KESULTS OF Sheer erie pry cOrMrNE 

Figures 9 and 10 show the vertical distribution of 
horizontal and vertical velocities which result from the 
numerical solution of eq. (54) for a given density distribu- 
meron i300, C= 2ttieeTncrtrmct £Oour normal modes are shown 
in each figure. Velocity values plotted are not absolute 
values but are relative to the horizontal velocity at the 
bottom for each mode. The ratio of w(z) to v(D) was computed 
using o and k values corresponding to a wavelength of 30 km. 
As can be seen on these graphs, the vertical velocity would 
Meme nsigniticant for a horezontal velocity on the order of 
cm/sec or larger. This means the only measurable motion 
weuwrd be in the horizontal dmrection. 

Figure 11 shows the effect of a change in the density 
@eceribution layer depth Gm the horizontal velocity distribu— 
tion. The second mode only is shown since effects on other 
modes are similar. The values chosen, 20 and 50 meters, 
represent the normal seasonal variation in cna interes of 
the under-ice layer. As can be seen, very little change in 
eaeevelocity distributionsresults. 

Figure 12 shows the effect on the second mode of a change 
in the curvature (e-fold depth, c) of the helow layer density 
distribution. This shows that as the e-fold depth increases 
(decrease in curvature), the velocity peak is less marked 
and occurs at a deeper depth. 

Frequencies were solved for hy using eq. (53) with an 


assumed horizontal wavelength of 30 km. Results were on the 
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Figure 9. Vertical Structure of Theoretical Horizontal 
Velocity Patterns for Modes 1 through 4. 
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Figure 10. 
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Srder of toes e665 or less, which corresponds to periods of 


tens of years. These values indicate the current patterns 

are essentially stationary. As can he seen from eq. (53) 
either a decrease in wavelength (increase in k) or an increase 
in mode (higher A) will cause the frequency to decrease. 

Thus the limiting, highest possible frequency for a given 
density distribution and horizontal wavenumher is for i = 0, 


a barotropic Rossby wave mode, and is given by: 


a S. 
on 


For a wavelength of 30 km this corresponds to a period of 
aaout 1O years. 

The horizontal distribution of velocities was arrived at 
Py assuming a length of 15 km to the first zero of the 
Semiicx-Ot) function of eq. (50) ane forming the product of 
ens trigonometric function and the 2Z(z) function at a con- 


stant depth. This is shown in Figure 13. 


eG. COMPARISON Ui diele OBSERVED DATA 

A comparison of the theoretical normal modes plotted in 
Figure 9 with the observed data in Figure 2 shows that the 
second normal mode most nearly matches the measured vertical 
M@eEstributcivon Of horizontal velocity. The only real discre- 
pancy between the two is the relative maximum in the theo-~ 
retical distribution at the surface. This appears to be an 
inherent characteristic of the CReore econ solution and is 
seen in each of the normal modes. One possible explanation 


is that these theoretical maxima were caused hy treatment of 
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Prource 1S. torimental Distmputionvwe: Theoretical 
Horizontal Velocities. 
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the under-ice surface as flat and frictionless in the develop- 
ment. The ice actually has a rough surface capable of dis- 
Sipating some energy of relatively slow moving currents. 

A second explanation is that more than one mode exists 
with a progressively decreasing amount of energy in succes- 
sively higher modes. Thus if Vv, (4) represents the horizontal 
welocity distribution for a given mode, then: 

eo 
v(z) = om av, (4): 


where a Pa 
n n+l’ 


would describe the actual velocity distri- 
bMtIOn. Assuming further, as a simplistic first approach, 
finait the only significant a, are ferme n= and n=2 (1,62, only 
the first two modes contain significant energy) the linear 
combination of modes 1 and 2 was plotted in Figure 14. The 
ay and a, were sized to give zero velocity at the surface, 
meaning a, was less than ay: Figure 14 can he seen to he 
closer to the observed data than the second mode alone. 
Obviously a more rigorous solution for the coefficients could 
be used to get an even closer fit to the data. This, however 
would not be physically meaningful. The purpose of this 
initial investigation was merely to show that the theoretical 
development would allow motion similar to what has been 
observed. This has been shown. The next phase of the inves- 
tigation, which will be discussed in Section IV of this 
thesis, is to study the initial value problem of how this 
specific ocean system responds to the initial perturbations 


in its density structure caused hy a newly forming lead 


(i.e., which normal modes, if any, are excited). 
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yo MM@ Neale (peweic(c) “elojijot sais) esta ll “euns, ie talciereeanela lh toul sy ciaalesvieafeinis 
more closely, a vertical cross-section through the origin 
was drawn showing theoretical horizontal velocity values, 
Figure 15. This was done for the second normal mode for a 
density distribution having h=30 and c=120, and for the 
PeMipne worl Zontal distributiom Gf Figure 13. Thus Figure 15 
could be compared with Figure 4, which is drawn for real 
eurrent data having a Similar density distribution. A 
Similarity can be seen in the general flow direction and 
in the shape of the isolines. Apart from the relative 
smoothness of the theoretical plot, the only Significant 
difference is in the upper homogeneous layer, which has heen 
discussed previously. Figure 16 is the same cross-section 
drawn for the summation of the first two normal modes. As 
can be seen, the cell structure is closer to the observed 


data of Figure 6 than was the second mode alone. 
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tee CE NERAL METHOD 

The problem addressed here is whether the density field 
perturbations which occur at a newly forming Arctic lead can 
force the theoretically developed oscillations to occur and, 
if so, whether the forced oscillations are consistent with 
observed data. 

The solution to this problem lies in formulating expli- 
citly the initial effect of the lead-induced perturhations 
on the pressure field. This effect can then be used as an 
MaereLal conditvon on the pressure function, ¢, in order to 
fully solve the Sturm-Liouville prohlem developed in Sections 
If and III. It was shown that the general solution to that 
problem was a set of eigenfunctions which gave the initial 
fjecponse of the Arctic Ocean to a haroclinic perturbation. 


Mie linear Combination of these functions can be written as: 


b (xe ee A Ban (2) COS (Kx 9 LE) (56) 
n=] m=l 
pemtchnat at Eime Zero: 
(812) a ALB alin (2) COs (Kk *) (57) 
ei 


MePmust be Moted that the SUMMation Of COlULmOnNs LNndicaced =oy 
eq. (56) is not necessarily a solution itself. It would only 


constitute a solution if the initial equations were actually 


Syl 





linear, which they are not. It appears reasonable to assume 
however that the non-linear contribution is small compared 
to the linear so that this summation approximates an actual 
solution. A method will now be shown whereby the values of 
AB Kn and Omn can bey LOuUncmmasea on, the initial specifica 


rom Of the ve field. 


IB. SeECreiCA Tre Cb ei EAE CONDITIONS 

The perturbations at a lead during formation and refreez- 
ing were described as small downward velocities and an 
increase in density under the lead due to salt rejection and 
small downward velocities due to shear processes. Consider- 
ing the increase in density below the lead as the dominant 
mreocess, the initial density perturbation will result in a 
deviation from the hasic state as shown in Figure 17; where 
4 shown is the depth of penetration of the density pertur- 
bation. Assuming that this takes place only directly under 
jae lead, thenma vertical cross-section Of initial density 
would appear as shown in Figures’. Using the hydrostatic 


assumption, this will result in an initial overpressure con- 


fbeton (perturbation pressure Po = 6) feechown in Figure 19 
fOr = os Be and zero elsewhere. This can then he used 
fepan initial Condition LO whicheeq. (5/)ecan be fitted in 


order to solve for the appropriate wavenumbers, frequencies, 

ema COnStants. A Lurther assumption that Z| <r ere icnbalbele 
h 

be made, thus a #"“=h can he used in describing the initial 


CrSwri pit von analytical iy. O04 (X12) Gan thus be written: 
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Figure dy. Pertuseariom Density Due to Salt 
Rejection at Time Zero. 


Figure 18. Vertical Cross-section of Perturbation 
Density at Time Zero. 
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Ix[ < x 


Z p 
dea = — = (58) 
od n OC Fae, 
= cee. 
barcoe. 
= e ip 
Z Pe 4 
p 
= 0 


([x] > x,) 


Ce. SOLUTION OF THE PROBLEM 
mine Initvaieeecondition, o0 (%12), 1s separable and can be 


written as 5 (x12) = f£(x)g(z) where: 


Ex) 0 els ee x) : (59) 


ll 
\—! 


Ix} < x) 


2710 (252 =i) (60) 


Q 
N 
T 


IL (27a) 


Thus, coefficients can be solved for separately from: 


t (ues ae BCOSs (kx) 
(61) 
g(z) = tant AZ, (2) 


Considering Ba first, these are the coefficients for a 


Fourier coSine series, which on the interval (0,L) are given 


iDy : 
x 
L x : 
weer Ee i firs wo) 
LB aa If £(x)dx + a IVE) ce 3 (62) 
0 x 
Ip - 
= = é ~ mS 5) We = 
Blo = f £ (x) cos | 4 | ax (nS ee) 
a eve <<’) iM=ae, Ces ) (63) 
mi m p ti fd ® @ 6 


where Kn = 'il/ tel sowerdtstance sutivetenitiy far removed) tron 
the lead so that ¢6 can be assumed to have a zero value. 

The values of B correspond to the zeroth mode, which is 
a barotropic mode that may exist but is not of interest in 
mime present Giscussion. PFPieure 20 is a plot of the approx- 
imation to £(x) formed by evaluating eq. (61) for m= ] 
through m = 40. It can be seen to approximate f(x) fairly 
elosely. Figure 21 1s a plot of the horizontal distrihution 
of v for these same modes. 

For each I thus found, and each eigenvalue corresponding 
to the Z (Zz) Vette aot lett mi nC tonic each wecam be “found 
a value of Orn’ 


The values of AL are given by: 


=) 
if Gig)Z (zjyaz 
0 n 
A = SE (ae Oe.) (64) 
n = 5 
i [Z (2) ) maz 
0 n 


for a general seriesS expansion, assuming the Z (2) are 
Seenogonal functions with respect to a weighting function 
of one. The denominator can be normalized for all modes by 


mmaaing a Lactor Ch Sue meet mack 


f (ep 7 (zc) |) dzeee) 
0 ‘oumee) 
or: 
: —il, 
= D 2 
Ca ls (Z_(z)) ar| (65) 


as = becomes: 
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Figure 20. Fourier Series Fxpansion of Initial 
Wevizontla ler reccure HilGctrihucron, 
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Figure 21. Fourier Series Fxpansion of Horizontal 
Nee) Gy inS eiiee On ees Cl ieee canon. 





a g(2)C 2, (2) dz (Waa 2 ts ) (66) 
0 


for the normalized distributions. The computation of A, can 
be further simplified by writing g(z) as the sum of an 
average and a variable part. Since D> h, the average part 


is essentially: 


jhe only variable part is then: 


Soap 2) = il as 77 (z > —-h) (67) 


Thus, eq. (64) can he written: 


= 10, = at 


AL, = / C2, (2) d2 + t (l+2/h)C 2, (2) dz (68) 


bate, che first integral in eq. (68) is equal to zero hased on 
mies Original continuity relationship wged in deriving 4(z). 
From Figure 9 it can be seen that each Z (2) is constant for 


@epths above -h. Thus eq. (68) finally becomes: 


A, == h/2° [aC] (ilo ee (69) 


ieee values were solved for using the IBM/360 computer and 
were then used to form the series summation given by eq. (61) 
for the first six modes. The results are shown in Figure 22. 
The approximation is not as close as the horizontal case due 
to using only six modes. This was a limit imposed by the 
numerical method. 


A complete set of values for Awe Beepeice, atic an is then 


mM mM 


available. Values for the distribution having c=120, h=30, 
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are given in Tables I, II, and III. Values of L=20 km and 
it km were assumed. 

The solution for $¢? given by eq. (56) will then begin 
Wet a Gistrimugeion as given by eq. (57), which is the ini- 
tial pressure distribution. Due to the extremely small 
values of o there will be no significant movement with time 
of the overall wave system. Phase velocities (o/k) are on 
fjne order of a Gmy7Sec. Pie Gictmdourion OF herizonte | 


velocities however, is proportional to: 


Cc 


VAIO 2 Se > AB men (2) Sink x-o, tt) (70) 
n=1 m=l1 ; 


mtewinitial horizontal and vertical distributions of v(x,z,t) 


were shown in Figure 21 and Figure 22. 
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Vv oeecONCLUSIONS 


The stated objective of this study was to model certain 
observed Arctic Ocean internal current patterns through an 
analytic study of the baroclinic response of the ocean to 
lead formation and refreezing. The problem was then divided 
mico two Major piasces: (l)™iinding the baroclinic internal 
normal modes of the Arctic Ocean; and (2) investigation of 
the forcing of these modes by the effects of lead formation. 

The normal modes of response found as a result of the 
theory developed in Section II showed a velocity pattern 
which was essentially stationary with no significant vertical 
velocity. However, horizontal velocities of the magnitude 
measured in the Arctic could exist within these stationary 
patterns. It was noted that the steepness of the pycnocline 
(measured by the parameter, c) had a direct and significant 
effect on the depth of maximum velocity, whereas the effect 
of a variation in layer depth was relatively minor. The 
apparent consistency of the cell structure seen in Figure 16, 
mime theoretical vertical cross-section of horizontal veloci- 
ties, with the observed data lends strength to the plaus- 
ibility of the developed theory. 

The fact that a simple linear combination of the first 
two theoretical normal modes, Figure 14, qualitatively 
Macches Meadcittecd velocity profiles a5 Of particular note. 


Pci CVemenemmeySles used to develop the theory to be 


G2 





consistent with actualsdata, in that obsexved current 
patterns appear to be part of the set of allowed solutions. 
Tt has not yet been shown conclusively that the perturbations 
at a lead will force these modes. 

The study done here is a first step which looked at the 
dynamics of a single Arctic lead. The theory developed has 
been shown to be dynamically possible hut the study has not 
yet looked at the energy involved in the system. It is 
possible that an interaction between leads or between lead(s) 
and the mean flow is needed to provide the energy required 
for the observed current magnitudes. 

The very small values of frequency, corresponding to 
periods of greater than 20 years, are attributahle to the 
magnitude of 8 in the Arctic. The vorticity equation derived, 
eq. (38), is a non-linear time dependent parabolic equation. 
However, the results found in the initial value problem 
studied in Seetion IV of this thesis indicate that the time 
dependence of the equation is not significant since the solu- 
tions were essentially stationary. Under the same geometry 
the problem could be reformulated as an emiipinc, Lincar 
equation forced with time dependent boundary conditions. 

An advantage of this type formulation would be the ability 
to write a separable differential equation in cylindrical 
coordinates, which was not possible using the development of 
Section II. Motion approximating the vortex-like patterns 
observed by Coachman and Newton [1973] could then he more 


easily described. 
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The plane wave solution for velocities developed in this 
thesis is a result of the formulation of the problem in 
rectangular coordinates with an infinite length lead. This 
worked well in developing the theory and resulted in a real- 
istic vertical section of velocity structure. However, 
actual leads are of finite length which could be expected to 


induce curvature in the velocity patterns. 
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VI. KECOMMENDATYONS "POR FURTHER RESEARCH 


The consistency found with observed data warrants the 
continuation of this study using the theory developed in this 
thesis as a hasis. The following recommendations for 
further research are made: 

im the feemalation of the initval value study as an 
elliptic problem, as discussed in section V; 

2. the use of vertical stretching caused by bathymetry, 
mie fLOrcea Convection at a lead, rather than the 6 effect as 
the driving mechanism in the developed vorticity equation; 

3. the investigation of the question of whether these 
velocity patterns have any influence on the stahility of the 
larger scale mean flow, including the effects of baroclinic 
instabilities; 


4. a study of the energetics of the system. 


65 





BIBLIOGRAPHY 


Seachi*an, bs K., Production of Supercooled Water During Sea 
ec VOBMaL tn smerec own arCtic Heat Budget and Atmos— 
Saeotlepe maculae Ines Rand Corporation, RM=5233-NSF, 
Mi AS Pe MSS 


Coachman, L. K. and J. L. Newton, Personal Communication, 
M97 3: 


Meerter, T. Dez, -Haline Convection Induced by the Freezing of 


Sea Wakeman Co OUNMSrEnheS even 757.5) 935-1908, March 
EG 8. 


mecter, T, De, sHaline Convection in Polynyas and Leads,” 
JOULE sree Geom. , V. 2, 402-469 ,, OGrtoper 1972. 


feeb ul. Aa, Current Measurements in the Canadian Basin of 
Phe Arctic. Ocean, University of Washington Department of 
Oceanography Technical Report 184, 1967. 


Goet, J. A., Ao Mamerical Investigation of Arctic Ocean 
Domamihess, wour. Phys. OGean., dNeplint, October 
1973a. 


fle, J. A., Unpublished Manuscript, 1973b. 


feelaebrand, F. B., Advanced Calculus for Applications, 
Prentice-Hall, Inc., 1962. 


foo, H. L., “Dynamics of Quasigeostrophic Flows and Insta- 


bility Gheery, Advances in Applied Mechanics, v. 13, 
2428-330 pereacemic Press, 1973. 


Peeelsen, Kaj ba, Methods in Numerical Analysis, MacMillan, 


1964. 

Bealosky, J.7eetne Stability of Currents in the Atmosphere 
PiCwocetwm——raite .,: sWOlUlrseAtmosc-. soci) vy. 2), 201-219, 
March 1964. 


Pedlosky, J., “Geophysical Fluid Dynamics," Mathematical 


Problems in the Geophysical Sciences, 49-60, American 
Mathemacreal Society, .lo7): 


Pilaccelanes yl eeehiiinetl Oa l =—sceudlves Orumc@envection CUrrentsS ate 
the Ocean's Surface Caused by Fvaporation, Radiation same 


a gene 











Oe ee 


Atmospheric Cooling, Final Technical Report for ONR 


Contract NO0014-67-A~ 0242-003, 1970. 


66 





Reid, KR. -O see Gee Me mebhee Crk lcCal Structurerot Macs ein 
Equatorial Wind-driven Currents of a Baroclinic Ocean,” 
JOUTr. Matimew Res... sue -5ol2, November 1948. 


Schaus, RR. Ho) Agihemmog mamvc Model of a Central Aretic Open 
Lead, M.S. Thesis, Naval Postgraduate School, 1971. 


Semcner, A, J., Personal Communication, 1973. 


67 





Lhe 


iT TAL bist RIBUTION 


Defense Documentation Center 
Cameron Station 
AlexanGd terme wag iinial 27:3 14 


Oceanographer of the Navy 
HoLiman Buamiidame No. 2 

2461 FEisenhower Avenue 
Alexandr tage virginia 22314 


Office of Naval Research 
Code 480 D 
Avr lIMG One ssc emma 222) 7 


Naval Postgraduate School 
Code 58 
Monterey, California 93940 


oie aiava eC OC Oa LZ 
Naval Postgraduate School 
Monterey, California 93940 


Dr. R. E. Stevenson 
Serene clemlaison OLfiice 


Scripps Institution of Oceanography 


ihaewollaerCalatornia 92037 


Dr. Jerry Ae Galt 


NOAA Pacific Marine Fnv. Lab. 
University of Washington WB-10 


Seattle, Washington 98105 
CDR Ricwawo lL. ltkin, USN 
96 Bonaire Circle 
SPULEGENeeew York 10901 


RK. Ke MeGregor 


Dinecteuencectic Program, Code 4)> 


Office of Naval Research 
Meainieeny Virginia 222107 


AIDJEX 

Division of Marine Resources 
University of Washington 
Seattle, Washington 98105 


68 


LoS. 


No. Copies 





el. 


Zs 


Assoc. Professor R. I. Williams 
Department of Meteorology 

Naval Postgraduate School 
Monterey, California 93940 


Professor L. K. Coachman 
Department of Oceanography 
University of Washington, WB-10 
Seattle, Washington 98105 


Ga 





UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (iKhen Date Entored) 








a e ok C 9A bf READ INS THUCTIONS 
1. REPORT NUMBER 2. GOVT ACCESSICN NO] 3. RECIPIENT'S CATALOG NUMBER 





- TYPE OF REPORT & PERIOD COVERED 
Master's Thesis; 
september 1973 


- PERFORMING ORG. REPORT NUMBER 


4. TITLE (and Subtitte) 






ieee tic MODEL OF OBSERVED SMALL 
seu BAROCLINIC CURRENTS IN THE 
ARCTIC OCEAN 


7. AUTHOR(s) 


Pemenmara ivan Itkin; 











SR TG PE Re ESTP L AOL Oe CPL ATO OD 


. CONTRACT OR GRANT NUMBER(S) 








LCDR USN 








. PROGRAM ELEMENT, PROJECT, TASK 


‘9. PERFORMING ORGANIZATION NAME AND ADDRESS 
AREA & VIORK UNIT NUMBERS 











Naval Postgraduate School 
Monterey, California 93940 














. REPORT DATE 


September 1973 


- NUMBER OF PAGES 





V4. 





CONTROLLING OF FICE NAME ANDO ADORESS 







Maver osctaraqcuate School 
MOntGereyyecalifornia 93940 


14. 








MONITORING AGENCY NAME & AOORESS(If different from Controlling Office) ». SECURITY CLASS. (of thie report) 





UneGlessciiace 





DECL ASSIFICATION/ DOWNGRADING 
SCHEDULE 







16. DISTRIBUTION STATEMENT (ol thle Report) 


meproved fOr public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (af tho ebotract entered In Block 20, (f different from Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reverao elde if neceeeary end identify by block nuaiber) 


ereriec; CURRENTS; BAROCLINIC; QUASI-GFOSTROPHIC 


20. ABSTRACT (Continue on reverse elde tf neceseary end identify by block number) 





An analytic model is developed of baroclinic currents 
measured by several investigators under the Arctic Qcean 
ice cover. - A time-dependent, Boussinesg, quasi-geostrophic 
VOrticity equation is developed to describe the flow. 
separable normal mode solutions are found which show 
Stationary current patterns qualitatively consistent with 


wR PEK SAI SIE EE SDT ETS OT ELIOT EE POP 











Satin ie Said Gore. 


DD , ORM 1473 EDITION OF 1 NOV 65 IS OBSOLETE 
(Page 1) S/N 0102°014+ 6601 | UNG hASS TE Las 











SECURITY CLASSIFICATION OF THIS PAGE (when Dsta Latcred) 


70 





UNCBASS IFIED 


SECURITY CLASSIFICATION OF THIS PAGE(Hhen Data Entored) 


observed data. An initial value prohlem is then formulated 
using density perturbations at a newly forming lead as a 
mercing Lunction. 


Se ——iS Oh ETO FE ETE OP LI DS 26 OLD. 





DD Form. 1473 UCTS) 
Pdan 72 We eoNool! | FD 
S/N 0102-014-6601 SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 


eal 








12 yunt? 





Thesis 14.5936 
~175 Itkin 
Cot An analytic mode! of 
observed small scale 
baroclinic currents in 
the Arctic Ocean. 





